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Abstract 

We present a Monte Carlo study of Kinetic Growth Walk on square as well as triangular lattice 
to show that it is not equivalent to canonical Self Avoiding Walk. 
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Introduction 

Growth is a natural process by which a system, for example, a tree acquires a history- 
dependent (or a causat) structure. Since a growing tree is never ungrown, its growth is an 
irreversible process. A generic but informal description of its structural features evolves from 
repeated observations of a variety of structural forms that these systems acquire during their 
natural growth. Whether it is also a meaningful statistical description depends on whether 
the collection of such irreversibly grown structures is amenable to a statistical mechanical 
study. In the case of causal and homogeneous networks [1] , it has been shown that growing 
(non-equilibrium) networks can be assigned appropriate growth probabilities and studied 
within the framework of equilibrium statistical mechanics. 

A trivial special case of a causal network is a Kinetic Growth Walk (KGW) [3-8] which 
may be considered as an ordered collection of 'nodes' (linear chain of monomers) such that 
two of them are of 'degree' one (end monomers) while the others are of 'degree' two each. 
It generates a Self Avoiding Walk (SAW) [2] configuration, step-by-step, by choosing one of 
the allowed steps at random with equal probability. Since the number of allowed steps at 
the growing end of the walk is a variable that depends on the configuration grown already, 
such a step-by-step growth is locally biassed. Averages taken over a collection of KGWs, 
uncorrected for the cumulative bias acquired during its growth, clearly establish the fact 
that KGW belongs to the same universality class as SAW [3-8]. In other words, like in the 
case of causal networks, the distinction between equilibrium (SAW) and non-equilibrium 
(KGW) configurations does not affect the scahng behaviour of the walk. This prompts us 
to take a closer view of how a collection of KGW configurations might be related to an 
(equilibrium) ensemble of SAWs. 

In fact, the bias-uncorrectcd KGW on a honeycomb lattice has been shown [9] to be 
exactly equivalent to a canonical ensemble of Interacting SAWs (ISAWs) [10] at the bath 
temperature proportional to {ln2)~^. These facts, put together, could be misconstrued to 
imply a general equivalence of KGW with ISAW at some bath temperature. Except for a 
hint [11] that kinetic walks may not be equivalent to canonical walks, no definitive study 
of this equivalence relation is available in the literature. In this paper, we present Monte 
Carlo evidence to show that such an equivalence does not exist on a square as well as on 
a triangular lattice. This is an explicit demonstration of the fact that the KGW is not, in 
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general, equivalent to ISAW at any bath temperature. 



A. Definitions 

Starting from some arbitrarily chosen site on a regular lattice of coordination number 
there are z possibilities for making the first step. For all subsequent non-reversal steps, we 
have a maximum oi z — \ possibilities to choose from. If, at any stage, we choose a step that 
takes us to a site already visited, we discard the process and start afresh from the beginning. 
This ensures that every site in the walk is visited once, and once only. It is clear that the 
probability of generating a successful A'"-step walk configuration can be written as 

where Zfe (= 1, 2, 3, • • • , 2; — 1) is the number possibilities available for the A;*'* step. If Zk is 
taken to be z — 1 for all the steps, then all the walk configurations are gernerated with equal 
probability, say PsAwiJ^)-, and are said to form an ensemble of Self Avoiding Walks (SAWs) 
[2]. 

On the other hand, if Z}- is the actual number of possibilities available for the W'^ step, 
then different configurations of the walk, called Kinetic Growth Walk (KGW) [3-8], are 
generated with different probabilities. 
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Whenever z^ < z — 1, we encounter non-bonded nearest neighbors, called contacts; the 
number of such contacts encountered at the k*^ step is given by ruk = {z — 1) — Zk, which 
when added up give the total number of contacts, n, in the walk. 

B. KGW on a honey-comb lattice 

The probability of generating an A'"-step KGW configuration with n contacts is given by 

[9], 

Pkgw{N) = e^'^^'PsAwiN) (4) 



which follows from the fact that the walk has no choice but to take the only available step 
after making a contact. Since the total number of A^-step SAWs, 'i^sAwi^), can be written 
as 

N 

^SAw{N) = J29{n) (5) 

n=0 

where g{n) is the number of SAWs with n contacts, the total number of KGWs on a honey- 
comb lattice is given by 

N 

^KGWhc{N) = J2g{n)e-'^^-^^'-' (6) 

n=0 

Thus, the KGW corresponds to the canonical ensemble of ISAWs [10] realized at an inverse 
temperature, Phc — in2 if the magnitude of the energy per contact, e, is taken to be unity. 
At any other 'bath' temperature, /3, KGW is identical to an ISAW at a lower temperature 
, {P + Phc)~^- Algorithmically, we first choose an available site at random (KGW rule) and 
then assign the Boltzmann factor if a contact was made. 



C. KGW on a square lattice 

The number of sites available for the current step is cither one or two or three depending 
on whether, in the previous step, the walk had made a double contact or a single contact or 
no contact respectively. The probabihty of generating an A'"-step KGW configuration that 
has made rii single contacts and n2 double contacts during its growth is therefore given by 

Pkgw{N) = e'^i'"^3/')+"2'"(3)P5Aw(A^) (7) 

Note that contacts made by the A^*'* step, if any, is not accounted for by the growth proba- 
bility Pkgw{N). So, the total number of contacts made by the walk configuration is 



ni + 2n2 
(ni + 1) + 2^2 

711+2(712 + 1) 



no contact at last step 
single contact at last step (8) 
double contact at last step 



It is not necessary that all the configurations with the same total number of contacts will be 
grown with the same probability because, growth being a history-dependent process, they 
could realize different pairs of values {7ii, 712} for the same value of n. 
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For example, consider an 8-step configuration encoded as 00323212 [12] that has a total 
of three contacts {rii — l;n2 — 1). The same configuration when grown in the reverse 
direction, encoded 03010122, has only encountered single contacts {rii = 3; 77-2 = 0). 

The fact that the total number of SAW configurations with n contacts, g{n), may be 
written as a double sum, 

J](5([ni + 2n2]-n), (9) 

rai,n2 

suggests that the total number of KGW configurations with n contacts, h{n), may also be 
given a similar expression - namely, 

h{n) = J2 e-['*i(-^)^i+"^(-'^)^2l5([ni + 2712] - n) (10) 

ni,n2 

where Pi = ln{3/2) and P2 = (l/2)Zn(3). They may be interpreted as the inverse 'tempera- 
tures' associated with the single and the double contacts respectively if e (= 1) is taken to 

be the (additive) energy per contact. Since {P2 > Pi), a double contact is 'colder' than a 
single contact! 

Configurations that are identical but grown in opposite directions are also counted into 
the expression for h{n) by virtue of the fact that they correspond to different pairs {rii, 772}. 
It is therefore clear that we have a distribution of {rii, n2}'s (or equivalently, a distribution of 
dimensionless contact-energy, S{ni, 712) = niPi+2n2P2 arising due to the history-dependence 
of the growth process. 

A canonical relation between h{n) and g{n) can be defined only if an average inverse 
'temperature' P{n) for a configuration with n contacts can be defined: 

h(n) = e-"(-^)^(")5(7z) where p(n) = ^^^A + ^ ^ ^ P2 (11) 

n n 

where < 7ii > and < 712 > denote the average number of single and double contacts respec- 
tively for a configuration with n contacts and may, in general, depend on n. Averaging out 
the history-dependence this way leads to the following expression for the total number of 
KGW configurations on a square lattice: 

N 

^KGWs,{N) = Yl 9in)e-<-'^^^^^ (12) 

n=0 

It is still not in the canonical form, in contrast to that for a honey-comb lattice (Eq.(6)), 
because P{n) is an intrinsic geometric parameter that characterizes a configuration with a 
given number of contacts, and not the inverse canonical (or 'bath') temperature. 
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FIG. 1: Size dependence of the average f3z for KGW. The asymptotic value of f3z is of the order of 
0.448 for a square lattice (Sq), and is of the order of 0.267 for a triangular lattice (Tr). 
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FIG. 2: The distribution of /?2 values for 300-step KGW configurations on a triangular lattice as well 
as on a square lattice. It is a truncated distribution to which a Gaussian could be fitted. The rele- 
vant parameters are (i) peak position ~ 0.265(rr); 0.446(5'(/); (ii) width ~ 0.0155(rr); 0.0203(S'g). 

A further averaging of the /? (n) values leads to the canonical expression 

N 

^KGWs,{N) = J2 5(n)e-"(-^)^- (13) 

n=0 

Identifying Psq as the inverse 'bath' temperature is, however, not justifiable because a canon- 
ical samphng of SAWs with Psq may not lead to KGWs. Such an interpretation is unique to 
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FIG. 3: The probability distribution for the number of contacts per step, n/N, of a KGW for 
N = 300. A gaussian could be fitted to this distribution for 300-step KGWs on a square lat- 
tice (Sq) as well as on a triangular lattice (Tr). The relevant parameters are, (i) peak value 
~ 0.351(59); 0.866(rr); (ii) width - 0.099(S'g); 0.1766(rr) 

honey-comb lattice only. 

The Monte Carlo estimates of Psq, plotted in Fig.l for various walk-lengths, N, extrapo- 
late to the asymptotic value ~ 0.448. The distribution of Psq values for a given walk-length, 
N — 300 in Fig. 2, is a priori expected to be a truncated one because Psq is bounded between 
/?i and (32- However, a Gaussian (width ~ 0.0203, peak position ~ 0.446) could be fitted to 
the distribution. 

The question remains whether (3^^ or (32^ could be identified as the 'bath' temperature. 
The fact that KGW belongs to the same universahty class as SAW does not provide a 
reason to a choose between the two, though they are both above the expected ^-temperature 
(~ 1.54). Moreover, such an arbitrary choice also leads to non-additive contact energy. 



D. KGW on a Triangular lattice 

A KGW on a triangular lattice is, in general, characterized by four inverse 'temperatures' 
{Pi = ln{5/4:),p2 = (l/2)/n(5/3),/33 = (l/3)/n(5/2) and p4 = (l/4)/n(5)) associated respec- 
tively with a single, double, triple and quadruple contacts. It is interesting to note that (3^ 
is greater than the expected ^-point value (~ 0.375), whereas the others are less. 
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The total number of KGWs on a triangular lattice may be given an approximate canonical 
form, similar to the one for KGWs on a square lattice (Eq.l3): 

N 

^KGWtriN) = J2 9{n)e--^-'^^- (14) 

where the asymptotic value of Ptr (~ 0.267) could be obtained by extrapolating the data 
presented in Fig.l. The truncated distribution of Ptr values, plotted in Fig.2 for walk-length 
= 300, is again fitted to a Gaussian, which peaks at ~ 0.265 and whose width is ~ 0.0155. 



E. Summary 

In general, a KGW on a two dimensional lattice of coordination number z is characterized 
by a set of {z — 2) inverse 'temperatures', /3m = 1, 2, ■ ■ ■ , z — 2), defined by 

f3^ = -ln( V ^, ^ J ; m = 1,2,--- ,z-2 (15) 



m \z — [m + 1] 
whose weighted average gives the canonical value 

N z-2 



A = E E (^) (16) 



n=l m=l 



where < > denotes the average number of m-contacts made. Note that (3^ = for walks 
without contacts (m = 0). The dimensionality dependence comes through the realizable 
values of n^s. As pointed out earlier, it is more appropriate to interpret Pz as an intrinsic 
geometric parameter that characterizes the KGW than to interpret it as an inverse 'bath' 
temperature. However, it may serve as a simple scale-shift parameter for the inverse 'bath' 
temperature if all the KGW configurations are reweighted so as to correspond to j3z. The 
reason why KGW turns out to be equivalent to SAW could be that the entire truncated 
distribution of lies below the corresponding 6'-point value. 
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